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t^. , 

j^-^ i Abstract. We present the results of the timing and color analysis of more than two hundred RXTE/PCA obser- 

^-j- . vations of the bright black-hole transient GX 339-4 obtained during its 2002/2003 outburst. The color-intensity 

evolution of the system, coupled to the properties of its fast time variability, allow the identification of four sepa- 
, rate states. Depending on the state, strong noise is detected, together with a variety of quasi-periodic oscillations 

at frequencies from 0.2 to 8 Hz. We present a characterization of the timing parameters of these states and com- 
^5 ' pare them to what has been observed in other systems. These results, together with those obtained from energy 

spectra, point towards a common evolution of black-hole transients through their outbursts. 

O ' Key words, accretion: accretion disks - black hole physics - stars: oscillations - X-rays: binaries 

c3 . 

• • . 1. Introduction an optically-thick but geometrically-thin accretion disk 
^ ' (Shakura & Sunyaev 1976). Very little short-term vari- 
k> ; Since the early years of X-ray astronomy, black-hole can- abmty in thc form of a p OW er-law shaped noise compo- 
. didate X-ray binaries have been known to show transi- nent; ig observed (see e.g. Belloni et al. 1999). The sec- 
CG tions between different spectral states. When X-ray in- Qnd ' ig the Low/Hard State ( LS) . the energy spec trum is 
strumentation became sufficiently sophisticated to allow dominated by a hard component, often associated with a 
detailed timing studies on short time scales, the defini- Com ptonizing region in thc accretion flow, with a high- 
tions of states were refined to include fast timing proper- energy cutoff energy of _ 10Q keV _ The pQwer spectra 
ties (see van der Klis 1995,2005; McClintock & Remillard can be deC omposed in a small number of flat-top (or 
2005). The number and defining properties of these states broad) com p 0ne nts, approximated with Lorentzian func- 
have changed with time (see e.g. Homan et al. 2001), but tionS; plug at timeg a i ow _f reqU ency quasi-periodic oscilla- 
it is now clear that fast timing variations are a key in- tion (Q p 0; sce Nowak 2000; p saltiSj Belloni & van der Klis 
gredient which needs to be considered in order to have a 1999; Belloni, Psaltis, & van der Klis 2002; Pottschmidt 
complete view of the states and state transitions. et aL 2004) ^ The characteristic frequencies of these corn- 
Three states are relatively well identified (see p0 nents follow rather precise correlations. Low-luminosity 
McClintock & Remillard 2005; van der Klis 2005). In neutron star X-ray binaries show very similar power spec- 
the High/Soft State (HS), the energy spectrum is dom- tra (see van der Klis 1995; Olive et al. 1998; Wijnands & 
mated by a soft thermal component, modeled with a van der Klis 1999; Psaltis, Belloni & van der Klis 1999; 
disk-blackbody (Mitsuda et al. 1984) with a typical tern- Belloni, Psaltis & van der Klis 2002; van Straaten et al. 
perature of -1 keV. This component is attributed to 2002,2003), with frequencies following the same correla- 

tion, indicating a probable common origin. The third is 
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els and characterizes the long periods of quiescence of 
transient systems. The timing properties here are not well 
known, given the low flux, although observations of XTE 
J1650-500 at a luminosity level of ~ 10 34 erg/s still sug- 
gest the presence of Low/Hard State timing properties 
(Tomsick, Kalemci & Kaaret 2004). 

In addition to these states, other behavior is observed 
in transient systems, often for much shorter periods. Their 
properties are complex and a simple classification was not 
possible until now. It is important to note that these ad- 
ditional states arc not only associated to the transitions 
between LS and HS (sec Homan ct al. 2001; Belloni 2004; 
McClintock & Remillard 2004; Casella et al. 2004). These 
states feature the largest variety of spectral and timing 
features. In the 1990's, a state associated to the highest 
luminosity observed in two systems, GS 1124-683 and GX 
339-4 was called Very High State (VHS; Miyamoto et al. 
1991,1993,1994). In this state, strong low-frequency (1- 
10 Hz) QPOs were observed together with noise variabil- 
ity showing rapid transitions between band-limited and a 
strong power-law power-spectral shape, while the energy 
spectrum is composed of both a hard and a soft compo- 
nent, with varying relative contributions. Two different 
QPO behaviors were observed depending on the energy 
spectrum, with fast transitions between them (Takizawa 
et al. 1997; Miyamoto et al. 1991). Overall, two separate 
"flavors" of VHS were reported: a hard one, dominated 
by the hard component in the energy spectrum and show- 
ing band-limited noise in the power spectrum, and a soft 
one spectrally dominated by the soft component and with 
power-law noise in the power spectrum. 

Later, some VHS properties were also observed at 
much lower luminosities and a new state was proposed, 
called the Intermediate State (IMS: Mendez & van der 
Klis 1997; Belloni et al. 1997). Homan et al. (2001) showed 
that in XTE J 1550-564 VHS properties were found at sev- 
eral well-separated luminosity levels and they proposed to 
consider the VHS as just the highest instance of the IMS. 
Both HS and LS were found to occur over a wide range 
of luminosities, effectively removing the need to distin- 
guish between VHS and IMS; Homan et al. (2001) pro- 
posed to consider the VHS as just the highest luminosity 
instance of the IMS (see also van der Klis 2005) . Recently, 
McClintock & Remillard (2005) proposed a new scheme 
for states of black-hole candidates with a somewhat differ- 
ent naming convention: the VHS was replaced by a more 
narrowly-defined Steep Power-Law state (SPL) and by 
several types of Intermediate States based on the parame- 
ters of the fitted spectral models. As noted by McClintock 
& Remillard (2005), there is no one-to-one correspondence 
between these primarily X-ray spectrally defined states 
and the time variability states we use here. 

Clearly, although there is general agreement about the 
LS and the HS as to their definition and their basic prop- 
erties, the states that differ from those two are complex 
and more difficult to classify and to interpret. Notice that 
the LS and HS are rather stable states, which can be ob- 
served for months or years, while the remaining states are 



often associated to state transitions and show strong and 
fast variations in their properties (see e.g. Miyamoto et 
al. 1994; Takizawa et al. 1997; Homan et al. 2001; Nespoli 
et al. 2003; Casella et al. 2004). Also, the detection of 
high-frequency QPOs (>20 Hz), although few are known, 
seems to be limited to these states (Morgan, Remillard 
& Greiner 1997; Remillard et al. 1999; Cui et al. 2000; 
Homan et al. 2001; Miller et al. 2001; Strohmayer 2001a,b; 
Homan et al. 2003a; Klein- Wolt, Homan & van der Klis 
2004; Casella et al. 2004; Homan et al. 2005b). It is also in 
these states that narrow and often strong low-frequency 
(<20 Hz) QPOs are observed. 

These QPOs have been classified into three separate 
classes. Wijnands ct al. (1999) and Homan et al. (2001) 
reported on two different types of QPOs in the RXTE 
data of XTE J1550-564: a broad one (type A), with a 
quality factor Q (the QPO frequency divided by the QPO 
full- width-half-maximum, FWHM)) of less than 3, and a 
narrower one (type B), with a Q larger than 6. Both QPOs 
were characterized by a centroid frequency of 6 Hz and as- 
sociated with a weak red-noise component, but with differ- 
ent phase-lag behavior. XTE J1550-564 also showed the 
more common QPO-type associated with a flat-top noise 
component (see Cui et al. 1999 and Sobczak et al. 2000). 
Remillard et al. 2002 dubbed this QPO 'type C: its fea- 
tures are a high coherence (Q > 10), a variable centroid 
frequency (in the range 0.1 - 10 Hz) and the simultane- 
ous presence of a strong fiat-top noise component (~10- 
40% rms). The type-C QPOs are seen to follow the same 
global frequency correlations as those of the low/hard 
state (Psaltis, Belloni & van der Klis 1999, Belloni, Psaltis 
& van der Klis 2002). While type-C QPOs are observed 
in many systems, the other two are less common. In ad- 
dition to XTE J1550-564, type-B QPOs were also ob- 
served in GX 339-4 (Nespoli et al. 2003), GRS 1739-278 
(Wijnands et al. 2001), XTE J1859+226 (Casella et al. 
2004), H 1743-322 (Homan ct al. 2005b) and possibly in 
4U 1630-47 (Tomsickfe Kaaret 2000), while type- A QPOs 
were observed in GX 339-4 (Nespoli et al. 2003), XTE 
J1859+226 (Casella et al. 2004), H 1743-322 (Homan et 
al. 2003b,2005b), and possibly in 4U 1630-47 (Tomsick 
& Kaaret 2000; Dieters et al. 2000) Furthermore, in the 
light of the A-B-C classification, the two QPOs observed 
in GS 1124-683 (Takizawa et al. 1997) and GX 339-4 
(Miyamoto et al. 1991) can be tentatively identified with 
types B and C, although a detailed re-analysis of Ginga 
data is necessary to confirm this association. These oscilla- 
tions, whose nature is still not understood, could provide 
a direct way to explore the accretion flow around black 
holes (and neutron stars). In particular, their association 
with specific spectral states and the phenomenology that 
is emerging indicate that they are a key ingredient in un- 
derstanding the physical conditions that give origin to the 
different states. 

State transitions and non LS/HS states appear to be 
important for the understanding of accretion onto black 
holes. There is evidence that the ejection of relativistic 
jets takes place during some state transitions (see Fender, 
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Fig. 1. Top panel: PCU2 light curve of the full outburst of 
GX 339-4, with one point per observation. Bottom panel: 
corresponding hardness ratio (see text). The dotted lines 
indicate major state transitions (see text). The horizontal 
dashed line shows the hardness level corresponding to the 
transition from the right branch to the left branch in the 
HID (see text and Fig. 2). 



known, with a lower limit of ^6 kpc (see Hynes et al. 
2004). 

GX 339-4 was the first BHC to show a radio/X-ray 
correlation in the LS (Hannikainen et al. 1998; Corbel et 
al. 2003; Markoff et al. 2003). Radio observations during 
the 1999 HS showed clear evidence of a strong decrease of 
core radio emission during this state (Fender et al. 1999b). 
During the 2002/2003 outburst, near the transition to the 
VHS (see Smith et al. 2002c), a bright radio flare was 
observed (Fender et al. 2002), which led to the formation 
of a large-scale relativistic jet (Gallo et al. 2004). 

It is clear that GX 339-4 is a very important source 
for our understanding of the accretion and ejection prop- 
erties of stellar-mass black holes. In this paper we present 
and discuss the results of the timing and color analysis 
of 205 RXTE/PCA observations of GX339-4 during the 
2002/2003 outburst. Preliminary results were presented in 
Belloni (2004). The results of a complete spectral analysis 
is presented in a forthcoming paper, while the analysis of 
the optical/IR/X-ray correlations is presented in Homan 
et al. (2005a). Spectral analysis of these data up to 2004 
Jan 29 was included in the long-term study by Zdziarski 
et al. (2004). 



Belloni & Gallo 2004, Corbel et al. 2004) and it is during 
these intervals that structural changes take place in the 
accretion flow. 

GX 339-4 is a transient black-hole candidate that is 
known to spend long periods in outburst. Historically, it 
was found prevalently in a hard state, although several 
transitions were reported (Maejima et al. 1984; Ilovaisky 
et al. 1986; Miyamoto et al. 1991). The system was one of 
the two that showed all 'canonical' X-ray states of BHCs 
(LS, HS and VHS: see Miyamoto et al. 1991). From the 
launch of the Rossi X-Ray Timing Explorer (RXTE) until 
1999 it remained bright, mostly in the Low/Hard state but 
with a transition to a softer state (see Nowak, Wilms & 
Dove 1999; Wilms et al. 1999; Belloni et al. 1999; Nowak, 
Wilms & Dove 2002; Corongiu et al. 2003). In 1999, the 
source went into quiescence, where it was detected with 
BeppoSAX at low flux levels (Kong et al. 2000; Corbel et 
al. 2003). A new outburst started in 2002 (Smith et al. 
2002a,b; Nespoli et al. 2003; Belloni 2004) and ended in 
2003 (Buxton & Bailyn 2004). After roughly one year in 
quiescence, a new outburst started, although it has not 
reached the same high luminosity levels to date (Buxton 
et al. 2004; Smith et al. 2004; Belloni et al. 2004; Homan 
2004; Kuulkers et al. 2004; Israel et al. 2004). A relativis- 
tically broadened iron emission line has been detected in 
the X-ray spectrum of GX 339-4, indicating the presence 
of a non-zero angular momentum in the black hole (Miller 
et al. 2004a,b). 

Recently, a high mass function (5.8±0.5 M Q ) has been 
measured for the system (Hynes et al. 2003), indicating 
strong dynamical evidence for the black-hole nature of the 
compact object. The distance to GX 339-4 is not well 
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Fig. 2. Hardness-Intensity diagram corresponding to the 
points in Fig. ^ Each point corresponds to one ob- 
servation. Dotted lines correspond to those in Fig. ^ 
Observation #0 is indicated by an empty circle. The gray 
symbols mark the observations included in the groups B 
(stars), A (diamonds), Ul (squares), U2 (circles) and weak 
(triangles) (see section 6). The gray arrows indicate the 
general time evolution along the outburst 
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2. Data analysis 

We analyzed a large set of 205 RXTE/PCA observations 
covering the 2002/2003 outburst of GX 339-4, partly from 
our own project and partly from data in the public archive. 
Notice that there are detections before MJD 52359 (see 
Homan et al. 2005a), but the source is too weak to de- 
tect signal in the power spectra. The observation log is 
shown in Table 1 (the table is limited to the observations 
mentioned explicitly throughout the text. The full table 
is available on line). The total amount of good data was 
555904 seconds. Since the observations correspond to a few 
different RXTE proposals, the observation modes differ 
between observations. We produced background-corrected 
count rates for each observation in the PC A channel range 
8-49 (3.8-21.2 kev), as well as a hardness ratio defined as 
the ratio of counts in the channel range 15-24 (6.3-10.5 
keV) over those in the channel range 8-14 (3.8-6.3 keV). 
As the data span a long period of time, we checked the sta- 
bility of the PCA gain by producing hardness and count 
rate from two archival observations of the Crab, from 2002 
March 28 and 2003 April 23. Between these two dates, the 
Crab count rate was found to increase by 0.02%, and its 
hardness to decrease by 0.2%. Both values are too small 
to affect our results. 

For timing analysis, we used custom timing-analysis 
software under IDL and MATLAB. For each observation 
we produced power spectra from stretches 128 seconds 
long using two separate energy bands: PCA channel band 
0-35 (corresponding to 2-15.3 keV) for the main power 
spectrum, and PCA channel band 15-40 (6.3-17.3 keV) in 
order to look for high-frequency oscillations, which are 
usually more prominent in this high-energy band (see 
Homan et al. 2002). 

The power spectra in the high-energy band were also 
used for a consistency check on the decomposition of 
broad-band noise into Lorentzian components. In the fol- 
lowing, for power spectrum we mean the low-energy power 
spectrum. For observations with ID 40031 the main power 
spectra were produced in the 8-35 (3.8-15.3 keV) channel 
range for lack of data at lower energies. The chosen time 
resolution was l/1024s, corresponding to a Nyquist fre- 
quency of 512 Hz, for the low-energy data and 1/8192 
for the high-energy data. Besides examining the spectro- 
grams, i.e. the time-frequency images, for each observa- 
tion we averaged all power spectra together. The high- 
frequency (> 256 Hz) part of the average power spectrum 
(where no signal could be seen) was fitted with a constant 
to estimate the level of the Poissonian noise. This aver- 
age was subtracted from the power spectrum, which was 
then rebinned logarithmically and converted to fractional 
squared rms. The power spectra were then fitted with 
a combination of Lorentzians (see Nowak 2000; Belloni, 
Psaltis & van der Klis 2002) using XSPEC vll.3. In the 
following, by flat-top Lorentzians we mean components 
whose centroid frequency is much lower than their width, 
often consistent with zero. 



We also produced cross-spectra between the data in 
channel range 0-11 (2-5 keV) and those in the channel 
range 12-29 (5-12.2 keV) from stretches of 16-s length. 
We then calculated averaged cross-spectrum vectors for 
each observation, from which we derived phase-lag spectra 
(see Casella et al. 2004). Positive phase lags indicate that 
the hard curve lags the soft one. In order to extract the 
relevant phase-lag information for different components, 
we also accumulated the phase lags in specific frequency 
ranges: 1-64 Hz to characterize the full spectrum, 1-6 Hz 
for the Lh component (see below), and for the QPO in 
a range centered on the QPO centroid frequency with a 
width equal to the FWHM of the QPO peak itself (see 
Reig et al. 2000). Notice that the estimate of the QPO 
lags depends on the relative level of the continuum and 
the QPO powers in the frequency range considered. When 
the continuum power is significantly larger than the QPO 
power, the lags of the continuum may well dominate. 

3. The outburst evolution 

In this section, we describe the general evolution of the 
outburst, as can be followed from Fig. ^ an( l ID hi Fig- 
^ we show the full light curve of the outburst and the 
corresponding evolution of the hardness. The Hardness- 
Intensity diagram is shown in Fig. Looking at Figs. Q 
and [3 we can identify four separate sections of the out- 
burst, corresponding to the four sides of the square traced 
by the source in the HID. The part of the outburst cov- 
ered here starts at the lower right. The source increases 
its flux steadily for 36 days, with only a small softening 
of the spectrum indicated by the slight slant of the right 
side of the square. Then the source on the HID starts 
moving rapidly to the left, with the source changing spec- 
tral shape but not its count rate. Notice that a constant 
count rate and a softening of the spectrum indicate that 
the source flux must be decreasing. The transition is in- 
dicated in both figures with a dotted line. The precise 
position of this line (between MJD 52394.44 and MJD 
52398.66) is determined by the timing properties (see be- 
low), but also by the properties of the optical/IR/X-ray 
correlation, which shows a marked reversal corresponding 
to this date (Homan et al. 2005a). This horizontal branch 
lasts about 10 days, after which the left branch is reached 
(after MJD 52410.53). The second transition is also shown 
with dotted lines: the precise positions of the transitions 
are determined by changes in the timing behavior (see be- 
low) . There follows a long period on the left branch, which 
is followed with larger relative variations in hardness than 
on the right branch. Also, the count rate is not monoton- 
ically decreasing on this branch. After an initial motion 
downwards, the source moves up and reaches a count rate 
similar to that of the first peak, near the top of the right 
branch (see also Fig. [IJ, after which it starts decreasing 
again. Finally, towards the end of the outburst (after MJD 
52693.73), there is a monotonic decrease in count rate and 
resulting motion downwards until GX 339-4 enters the 
lower horizontal branch, which eventually curves and be- 
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Table 1. Log of the PCA observations (limited to the ones 
explicitly mentioned in the text. The headers mark the four 
source states (see text). The full table is available only in elec- 
tronic form at the CDS.) 



Obs 


Obs. ID 


MJD 


Exp. 


PCU2 


HR 


# 






(s) 


c/s 




Low-Hard State 





60705-01-54-00 


52359.64 


768 


9.0 


0.91 


1 


60705-01-55-00 


52366.60 


896 


153.2 


0.88 


2 


70109-01-02-00 


52367.76 


6784 


235.3 


0.87 


8 


70110-01-04-00 


52381.11 


1536 


868.4 


0.82 


17 


70109-01-04-00 


52387.56 


10880 


956.0 


0.82 


24 


70110-01-08-00 


52394.44 


1536 


1003.3 


0.81 


Hard Intermediate State 


25 


70110-01-09-00 


52398.66 


1024 


1027.3 


0.79 


27 


70108-03-01-00 


52400.85 


8192 


1023.4 


0.74 


28 


70110-01-10-00 


52402.49 


1024 


1003.0 


0.59 


29 


70109-04-01-00 


52405.58 


5888 


1136.7 


0.37 


30 


70109-04-01-01 


52405.71 


17408 


1134.3 


0.37 


32 


70110-01-11-00 


52406.70 


896 


1132.3 


0.35 


33 


70110-01-12-00 


52410.53 


1152 


1081.7 


0.28 


Soft Intermediate State 


34 


70109-01-07-00 


52411.60 


5376 


1102.5 


0.24 


35 


70110-01-13-00 


52412.07 


896 


1033.4 


0.22 


41 


70110-01-14-00 


52416.60 


768 


1011.0 


0.26 


42 


70110-01-15-00 


52419.24 


896 


940.7 


0.26 


43 


70108-03-02-00 


52419.43 


10368 


827.2 


0.22 


94 


70110-01-45-00 


52524.95 


1280 


1161.3 


0.23 


95 


70110-01-46-00 


52527.85 


1024 


1102.8 


0.22 


96 


70109-01-23-00 


52529.58 


1792 


1062.7 


0.22 


97 


70110-01-47-00 


52532.75 


1280 


1051.8 


0.23 


98 


70110-01-48-00 


52536.11 


640 


957.7 


0.23 


99 


70109-01-24-00 


52536.36 


2816 


945.1 


0.21 


103b 


70130-01-02-00 


52546.51 


9344 


769.9 


0.19 


109 


70110-01-55-00 


52558.63 


1024 


564.4 


0.15 


High/Soft State 


110 


70109-01-26-00 


52560.41 


3200 


550.7 


0.12 


148f 


70110-01-86-00 


52693.73 


896 


82.3 


0.19 


Hard Intermediate State 


149 


70109-01-37-00 


52694.92 


3712 


87.2 


0.22 


149g 


50117-01-03-00 


52706.84 


7808 


78.9 


0.36 


149h 


70110-01-89-00 


52707.92 


896 


72.7 


0.26 


153c 


50117-01-04-00 


52715.85 


17408 


52.1 


0.14 


154 


60705-01-57-00 


52716.70 


2560 


59.0 


0.16 


154b 


70110-01-92-00 


52717.50 


1024 


60.6 


0.18 


155 


60705-01-59-00 


52731.56 


3456 


77.7 


0.62 


157b 


70110-01-01-10 


52740.01 


896 


56.4 


0.79 


Low/Hard State 


158 


80116-02-02-00 


52741.70 


6016 


51.8 


0.79 


158b 


70128-02-03-00 


52742.24 


14464 


49.2 


0.80 


161 


60705-01-61-00 


52746.97 


2816 


33.3 


0.84 



comes once more vertical, ending its outburst on a position 
very close to its starting position. In order to identify these 
branches in terms of 'canonical' BHC states, we need to 
investigate the timing and spectral properties. A spectral 
analysis will be presented in a forthcoming paper. In the 
following, we describe the detailed evolution of the power 



spectrum of GX 339-4 from a purely phcnomcnological 
point of view. The relation between our results and the 
source states, as determined also from spectral parame- 
ters, will be examined in the discussion. 




0.01 0.1 1 10 100 



Frequency (Hz) 



Fig. 3. Four average power spectra from the right branch 
and four from the top branch. From top to bottom, 
they correspond to observations #2,17,24,25,27,28,30,33. 
The power spectra are plotted in the vP v representation 
(Belloni et al. 1997) and are shifted vertically for clarity. 
The increasing trend at high frequencies that can be seen 
in some power spectra is made of 3<r upper limits, shown 
as inverted triangles. 



4. The right branch 

The first significant power detected in the power spec- 
trum corresponds to observation #0. Here, because of the 
low count rate, a fit with one Lorentzian gives a reduced 
X 2 close to unity, but the fit is clearly not good, as con- 
siderable excess power is visible at low frequencies. A fit 
with two Lorentzians gives as characteristic frequencies 
0.021±0.7 Hz and 0.32±0.08 Hz, with a total integrated 
fractional rms of ^25%. As these frequencies are much 
lower than the ones observed in subsequent observations, 
we do not know how to classify them and do not include 
them in the Tables. 

Observation #1 takes place one week after Obs. #0 
and the count rate is a factor of 16 higher. The power spec- 
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trum of Obs. #2 is shown in Fig. (top spectrum). The 
shape of the power spectrum is very similar to that ob- 
served by RXTE in GX 339-4 during its typical Low/Hard 
State (see Belloni et al. 1997; Nowak, Wilms & Dove 
1999) and also similar to that of most BHCs in that state 
(McClintock & Remillard 2005). From Obs. #1 (similar to 
Obs. #2) to Obs. #24, the power spectrum has a similar 
shape, with its characteristic frequencies increasing with 
time and therefore with count rate. All power spectra can 
be fitted with up to four Lorentzians: one at low frequen- 
cies (ii), one at intermediate frequencies (L2), one at high 
frequencies (L3) and a QPO peak (Lq, observed only in 
a few cases). For the first two observations, an additional 
Lorentzian component Li is needed between L\ and Li- 

The characteristic frequencies (defined as v maxi see 
Belloni, Psaltis & van der Klis 2002) of these components 
can be seen in Table 2. Their time evolution, with the ex- 
ception of the Li component, present in only two spectra, 
can be seen in Fig. 01 together with the tight correlation 
between v\ and 1^2- 

The integrated fractional rms of the L\ component de- 
creases smoothly from about 30% to ~17%, while at the 
same time the L2 component goes from 28% to 14%. The 
fractional rms of the L3 component, when present, is al- 
ways 2-3%, and the few detections of Lq are at a level of 
4-7%. The quality factor Q is always less than unity for 
L\ and L2, ~1 for L3 and typically ~8 for Lq. 
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Fig. 4. First three panels: time evolution of the charac- 
teristic frequencies of the L\, L2 and L3 components for 
the top/right-branch observations at the beginning of the 
outburst. Fourth panel: v% vs. v\ correlation. In all panels, 
the dotted line separates the right-branch points from the 
top-branch points. 

In order to identify the broad components with those 
presented by Belloni, Psaltis & van der Klis (2002), we 
need more detections of the QPO component. Therefore, 
this identification will be shown at the end of the next 
section. In any case, all the power spectra corresponding 
to this branch of the HID have a shape that allows us to 
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Fig. 5. Phase lags for the right and top branches, (a) Time 
evolution of the lags of the 1-64 Hz continuum, (b) Time 
evolution of the QPO lags, (c) Correlation between the 
lags of the 1-6 Hz continuum and the characteristic fre- 
quency v\ of the L\ component, (d) Correlation between 
the lags of the QPO and its characteristic frequency. 

classify these observations as belonging to the low/hard 
state. 

The phase lags for the continuum, for the QPO, and 
between 1 and 6 Hz can be seen in Fig. El All lags are 
positive. For the QPO there are only three points, but the 
other two show a clear increase in the value of the hard 
lag from less than 0.05 radians at the start to almost 0.15 
radians at the top of the right branch. A representative 
phase-lag spectrum (from Obs. #8) is shown in Fig. 
Hard lags are also typical of the low/hard state (see e.g 
Pottschmidt et al. 2000). 

5. The top branch 

Starting from Obs. #25, the hardness ratio starts to de- 
crease more rapidly (see Fig. |2J an d a clear QPO with 
harmonic content appears in the power spectrum. This 
can be seen in the bottom four spectra from Fig. El which 
also shows that the characteristic frequencies continue to 
increase. The shape of the power spectrum is typical of 
the Very High State (VHS) in its hard version with band- 
limited noise (see Miyamoto et al. 1991,1993,1994; Homan 
& Belloni 2005). 

In order to fit these power spectra, additional compo- 
nents such as harmonic peaks and QPO shoulders (see e.g. 
Belloni et al. 1997) are needed. We decided not to include 
the shoulders in our fits in order to use a model consistent 
with that used for the previous observations. The only ad- 
ditional component is a broad Lorentzian (Lq) appearing 
at low frequencies, while the L3 component becomes un- 
detectable, likely because its characteristic frequency in- 
creases and leaves the frequency range where we detected 
significant power. The characteristic frequencies are re- 
ported in Table 2 and their evolution is shown in Fig. 0] It 
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Obs # 




V\ 




V2 




Vi 


Right branch 


1 




a ni o i a aact 

0.012 ±0.005 




O OO 1 A OA 

z.ZoiU.zU 




a oo i a r\a 

0.28±0.06 


2 




n AO 1 A A 1 

0.03 ±0.01 




O OO 1 A A C. 

z.ZZ±U.U5 




A O O 1 A A A 

0.38±0.04 


3 




0.06 ±0.01 




O CO 1 A AT 

2.52±0.07 






4 




A 1 A Al 

0.07 ±0.01 




O A 1 (A AT 

Z.41±U.U7 






5 




A A£? 1 A Al 

0.06 ±0.01 




O O A 1 A AO 

2.34±0.08 






6 




A AO 1 A Al 

0.08 ±0.01 




O O O 1 A 1 

Z.88±U.zl 






-r 

1 

8 




A 1 O 1 A A 1 

0.13 ±0.01 

n in 1 A A 1 

U.lz ±(J.U1 




O 1 1 A A C. 

3.21±0.05 

O AA 1 A AT 

3.09±0.07 


O A 1 A _l_ 1 OC 

z4.14il.oo 




9 




A 1 A 1 A A 1 

U.14 ±U.U1 




O C\f~' 1 A AA 

3.06±0.09 


O A O CT 1 O O O 

z4.85iz.bz 




i n 




nir 1 A A 1 

0.15 ±0.01 




O AO 1 A AO 

3.08±0.08 






11 




A 1 A 1 A A 1 

U.14 ±U.U1 




O OC 1 A A f. 

2.85±0.05 






12 




A 1 A 1 A A CT 

0.14 ±0.05 




O AA 1 A A A 

2.99±0.04 


01 AO 1 O OO 




13 




A 1 O 1 A A 1 

0.13 ±0.01 




O C\G LA 1 C 

3.06±0.16 






14 




0.15 ±0.03 




O 1 1 A AO 

3.21±0.03 


OA O O _l_ 1 C 1 

zy.oozbl.ol 




i - 

lo 




A i 7 _l_ A A 1 

0.17 ±0.01 




O OO 1 A A a 

3.33±0.05 


OA AC\ I A OO 

o4.4L)±4.zo 




1() 




A OA 1 A AO 

U.zU ±U.Uz 




O A A 1 A 1 O 

o.49±U.lz 


OA OO 1 O OA 




1 1 




n i7 -i_n m 


U. IDztU.Ul 


O.ZOztU.Uo 


OQ qo_|_n Q7 

zo. ooztu.y ( 




18 




0.20 ±0.01 


0.20±0.01 


3.54±0.08 


28.87±3.66 




19 




0.16 ±0.01 




3.33±0.08 


33.34±3.26 




20 




0.18 ±0.01 




3.06±0.11 






21 




0.23 ±0.01 


0.20±0.01 


3.71±0.07 


36.48±2.74 




22 




0.22 ±0.01 




3.39±0.10 


30.08±4.60 




23 




0.22 ±0.01 




3.81±0.14 


26.15±2.68 




24 




0.26 ±0.02 




3.72±0.08 


38.68±3.14 




25 




0.55 ±0.05 




4.88±0.16 


35.20±2.81 




Top branch 


26 




1.21 ±0.03 


1.26±0.01 


9.59±0.23 


39.31±1.64 




27 




1.23 ±0.02 


1.33±0.01 


9.63±0.31 


43.65±1.85 




28 




3.31 ±0.12 


2.08±0.01 




49.33±6.03 




29 


1.08±0.01 


6.46 ±0.26 


5.46±0.01 








30 


0.99±0.01 


5.93 ±0.13 


5.45±0.01 




96.93±14.78 




31 


0.94±0.02 


5.73 ±0.61 


5.34±0.02 








32 


0.99±0.01 


6.08 ±0.75 


5.82±0.02 








33 


0.99±0.01 




8.12±0.18 









is evident that the components detected in the top branch 
are consistent with being the extension to high frequency 
of those seen in the right branch. However, a gradual tran- 
sition can be seen in the top right panel of Fig. Despite 
the fact that the energy spectrum changes considerably, 
as indicated by the hardness ratio, the timing properties 
do not seem to show sharp changes, with the exception of 
the appearance of a clear QPO of type C (see Miyamoto 
et al. 1994, Remillard et al. 2002 and Casella et al. 2004). 
The characteristic frequencies increase, and in particular 
V\ and Vi show a fast time evolution to higher values. The 
fourth panel in Fig. 0] shows their correlation. Figure 
shows the evolution of v\ for the right and top branches 
(filled circles) as a function of count rate and hardness. 
The gradual transition between the two branches is evi- 
dent. 

We can now attempt an identification of the compo- 
nents present in our power spectra and those presented 
in Wijnands & van der Klis (1999), Psaltis, Belloni & 
van der Klis (1999) and Belloni, Psaltis & van der Klis 



(2002). It is natural to identify the Lq component with 
the low-frequency QPO in BHC. The characteristic fre- 
quency of the L\ component is always close to vq and 
on that basis can be identified with the Lh component 
from Belloni, Psaltis & van der Klis (2002). Following this 
scheme, Lq corresponds to Li,, L2 to Li and L3 to L u 
(see Belloni, Psaltis & van der Klis 2002 for the descrip- 
tion of these components). In order to check whether these 
identifications hold, we plot the three global correlations 
(WK from Wijnands & van der Klis 1999; and PBK from 
Psaltis, Belloni & van der Klis 1999) in Fig. together 
with the corresponding points from our observations. It is 
worth noting that the Lj, (L ) component is present here 
in only a few observations, while it constitutes the dom- 
inant component in other systems (see Wijnands & van 
der Klis 1999). This is also seen in many power spectra 
of GRS 1915+105 (see e.g. Morgan, Remillard & Greiner 
1997). 

Due to this, we can put only four points on the WK 
relation. Also shown in Fig. |H1 (top panel) are the V\ — vq 
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Fig. 6. Phase-lag spectra for three representative obser- 
vations from the right (Obs. #8), top (Obs. #29) and left 
(Obs. #34) branches. The vertical dashed lines mark the 
centroid frequency of the QPO. 




100 

PCU2 counts/s 




Fig. 7. Evolution of the v\ frequency as a function of count 
rate (top panel) and hardness (bottom panel) for the ob- 
servations of the right and top branches (filled circles) and 
of the bottom branch (empty circles). The X axes of both 
panels are the same as on the corresponding axis on Fig. 
□ 

pairs, which all have v\ vq. This second correlation was 
not present in the original Wijnands & van der Klis (1999) 
work, but was presented by Belloni, Psaltis & van der Klis 
(2002). 

For comparison, from the early-outburst observations 
of XTE J1550-564, one can see that the and com- 
ponents have a rather different energy spectrum, with 
appearing much more prominent at higher energies (see 
Fig. 4 in Cui et al. 1999). However, a check from the high- 
energy power spectrum from our observations shows that 
L does not appear at higher energies. Fig. |S] suggests that 
the identification v\ = vh and V2 = vi is correct: in par- 
ticular, the alternative interpretation v\ — Vb and = vy L 
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Fig. 8. Top panel: correlation between v\ and vq (large 
circles) and points from the WK correlation (Wijnands 
& van der Klis 1999, Belloni, Psaltis & van der Klis 
2002). The upper points are our detections of v$ from 
the top branch. Middle panel: correlation between vq and 
V2 (large circles) and points from the PBK correlation 
(Psaltis, Belloni & van der Klis 1999, Belloni, Psaltis & 
van der Klis 2002). Bottom panel: correlation between 
and V2 (large circles) and points from the PBK second 
correlation (Psaltis, Belloni & van der Klis 1999, Belloni, 
Psaltis & van der Klis 2002). 



provides no match at all to the WK relation in the top 
panel of Fig. |H1 The agreement visible in Fig. |H1 shows 
that the identification is correct. 
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It is interesting to notice that the integrated fractional 
rms of the L\ component continues its smooth decrease 
observed in the LS, reaching ^8% at Obs. #32, while the 
rms of Li drops to a few % before the component becomes 
undetectable. The QPO has a rather stable rms between 
4-5%, similar to that of the few detections along the right 
branch. 

The evolution of the phase lags (see previous section) 
along the top branch, after the transition, can be seen in 
Fig. while one representative phase-lag spectrum (for 
Obs. #29) is shown in Fig. |SJ Once again, hard lags are 
observed. Also the lags seem to show a smooth evolu- 
tion through the state-transition. These phase lags are in 
marked contrast with what is usually observed in VHS 
observations, where the type-C QPO shows soft lags (see 
Wijnands et al. 1999; Reig et al. 2000; Homan et al. 2001; 
Remillard et al. 2002; Casella et al. 2004). 

The last observation of this state, Obs. #33, four days 
after Obs. #32, shows a higher QPO frequency (~7.8 Hz) 
which can be classified as type C* (see Casella et al. 2004). 
Only the Lq component can be detected here. An addi- 
tional band-limited component appears at lower frequen- 
cies, with a characteristic frequency ~0.16 Hz. 

In order to compare type-C QPOs with type-B QPOs, 
in Fig. El we show the spectrogram for Obs #30 (see Sect. 
6.1). No clear time variability of the centroid frequency of 
the type-C QPO is seen. 




v y> tK \m iw at in in *w tt an 

ftl'tEU'VjIiJjl 



Fig. 9. Top panel: a 500-s segment of the light curve of 
Obs. #30 (1 second bin size); count rate is total count rate 
for all PCUs. Bottom panel: corresponding spectrogram 
(time integration 4 seconds, time step 1 second), where 
darker points correspond to higher power. Inset: average 
power spectrum from the full observation. 



6. The left branch (top section) 

After Obs. #33, the power spectrum changes radically and 
no more type-C QPOs are observed until the source leaves 
the left branch on the HID. As one can see from Fig. 



Table 3. Log of the five types of power spectrum in the left 
branch. Grouped observations are indicated by '-' symbols 



Type 


Observations 


No signal 


65,68,72,73,76,77-78,83-84,89-90,100, 




105-106 


Red noise 


56,57-58,60,63-64,70-71,74-75,86 


Weak QPO 


44-46,66-67,69,79,82,85,88b,91-92,94,96, 




98-99,101-102,103b,107 


LF noise 


35,36-40,47-50,51,52-55,59,61-62,80-81, 




87,88,93,95,103,104,108,109 


type-B QPO 


34,41,42,43,97 



121 after having crossed the second dotted line, GX 339-4 
spends a long time to its left. During this period, the flux 
does not change monotonically (see Fig. ^ . This complex 
interval continues up to observation #109 (MJD 52558). 
After this time, the PCU2 count rate remains below 500 
cts/s and only a weak power-law signal is observed in the 
power spectrum: this bottom section of the left branch will 
be described in the next section. During the period before 
Obs. #109, the timing behavior of GX 339-4 is complex. 
We summed the power spectra of consecutive observations 
which were close in time and had a similar position in the 
HID. The resulting power spectra can be summarized in 
a few types (see Table 3): 

— no signal: no significant variability is detected in the 
power spectrum. 

— red noise: a featureless red-noise component can be 
seen. 

— weak QPO: typically, in the low-energy band, a broad 
QPO is detected at a frequency of ~5 Hz, together 
with a broad-band noise component with a character- 
istic frequency of a few Hz. In the high-energy band a 
QPO at ~10 Hz is also seen. All these features have 
an integrated fractional rms of a few %. 

— low-frequency noise: some low-frequency power in the 
form of a weak band-limited noise is observed, with no 
significant QPO detection either at low and high ener- 
gies. Sometimes a weak QPO is also detected, making 
the distinction between this type and the previous one 
uncertain. 

— type-B QPO: a narrow transient type-B QPO (sec 
Nespoli et al. 2003) is detected during the full obser- 
vation or only during a part of it. This case will be 
discussed in detail below. 

A pattern in this complex behavior can be found by 
integrating the fractional rms in the 0.1-64 Hz range for 
all power spectra and plotting the result as a function of 
hardness. The result, which includes also the points from 
the right and top branches, can be seen in Fig. ITol 

Here, right branch, top branch (type-C QPO) and 
type-B power spectra are clearly separated. Type-C points 
follow a correlation that extends all the way to the right 
branch points, while type-B points also follow a correlation 
that deviates from the type-C one, but possibly connects 
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Fig. 10. Plot of integrated fractional rms versus hard- 
ness for the observations of the right, top, and upper 
left branch. Different symbols correspond to: right branch 
(open circles), top branch (asterisks), type-B QPO (di- 
amonds), weak QPO (filled circles), LF noise (squares), 
red noise (plus signs) and no signal (dots). For the the 
definitions of the different elliptical regions, see Sect. 6. 



The full description of the results for these observations 
can be found in Appendix A. 















.""'yv<* ,"" 


ttt ) , t t '«tt , t "tt"ttV 


," ( y 






* i 






'V 




'tjttttf.tV'jtfMftWt 




t| tt r 




















' ' ' ' » ' f ttttj 




f 







Frequency (Hz) 



Fig. 11. The Poissonian-subtracted power spectrum of 
the five observations with a type-B QPO (see text), 
vertically shifted for clarity. From top to bottom: Obs. 
#34,41,42,43,97. 



to it. A few observations occupy a specific region of the 
diagram, roughly on the extension of the type-B corre- 
lation. In analogy with what is seen in XTE J1859+226 
(Casella et al. 2004) we call this region 'type A' (see be- 
low). The other points are also grouped: we identify a 
group including almost all observations with no signal, and 
two additional groups that we call 'Ul' and 'U2'. These 
three groups correspond to different levels of integrated 
rms (see Fig. 1101) . Ul and U2 are ketp as separate sub- 
groups as the power spectra of most of the observations in 
Ul are classified as "LF noise", while those in U2 show a 
weak QPO. This is reflected in the average power spectra 
(see below) . Notice that the observations corresponding to 
these groups are marked with different symbols in Fig. ■ 
In the following, we will examine the average power spec- 
trum of these groups, and in more detail the instances of 
type-B QPO detections. 

6.1. Type-B QPO 

The first detection of a type-B QPO in this outburst of GX 
339-4 has already been reported by Nespoli et al. (2003). 
The QPO in that observation (Obs #34) is transient: dur- 
ing the first part of the observation only a broad bump is 
present in the power spectrum, while the narrow QPO 
appears suddenly and shows marked variations in its cen- 
troid on a characteristic time scale of ~10 s (Nespoli et al. 
2003). A similar QPO was present also in XTE J1859+226 
(Casella et al. 2004). We investigated the time evolution 
of the QPO of the other four type-B observations by pro- 
ducing spectrograms with a time resolution of 4 seconds 
(see Nespoli et al. 2003) and found that indeed in all four 
cases the QPO is not present for the whole observation. 



In summary, (transient) type-B QPOs are found in 
five observations (see Fig. 111(1 : their centroid frequency 
ranges from 4.5 to 7.0 Hz and their fractional rms is in 
the range 5-8%. As in the case of the type-B QPO re- 
ported by Nespoli et al. (2003), the hardness increases 
slightly (by ^0.02) when the QPO is present compared to 
the no-QPO intervals. Selecting QPO- and no-QPO inter- 
vals (see Appendix A), we can also notice that for obser- 
vations #34,43,97, the no-QPO 0.1-64 Hz rms is around 
3%, while when the QPO is present it is around 8-9%. For 
the remaining two observations, #41,42, the no-QPO and 
QPO fractional rms is always around 9%, but as these are 
the observations where the QPO appears intermittently, 
our selection of the no-QPO interval was approximate and 
it cannot be excluded that some residual spurious power 
from the QPO intervals and/or from the rise and decay of 
the flux were included in the power spectra. 

An example of the phase-lag spectrum for an obser- 
vation showing a type-B QPO (Obs. #34) can be seen 
in Fig. H3 The errors are large, but around the QPO fre- 
quency clear hard lags are observed, in line with what was 
observed in previously reported cases of type-B QPOs (see 
e.g. Casella et al. 2004). 

6.2. Type- A QPO 

The average power spectrum of the observations in the 
'type A' group (Obs. #35,94,95,96,98-99,103b) can be 
seen in Fig. It can be fitted with a model consisting 
of two flat-top Lorentzians with characteristic frequencies 
0.07±0.02 Hz and 0.54±0.13 Hz respectively, plus a broad 
QPO with centroid 7.49±0.18 Hz and FWHM 3.64±0.56 
Hz. The integrated fractional rms of the three components 
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is 2.8%, 4.7% and 3.7% respectively. In the high-energy 
band, the power spectrum is qualitatively similar. The 
broad shape could be related to the averaging of differ- 
ent observations, but this shape is remarkably similar to 
that observed by Nespoli et al. (2003) in the interval be- 
fore the appearance of the type-B QPO. The similarity in 
shape, centroid frequency and total rms make us identify 
this QPO as of type A (see Wijnands et al. 1999; Casella 
et al. 2004 and references therein). 




0.01 0.1 1 10 

Frequency (Hz) 



Fig. 12. Average power spectra of the observations clas- 
sified as (from top to bottom) type A, weak signal, U2 
and Ul (see Fig. llO|) . The top power spectrum is in power 
units, all the others are shifted downwards by a factor 
of 10 from each other. The lines are the best fit models 
described in the text. 



6.3. Weak signal power spectrum 

The average power spectrum of the observations in the 
'weak signal' group is shown in Fig. ^| Its total 0.1- 
64 Hz fractional rms is around 1.3% and can be fitted 
with a power law (index 1.16±0.03) plus a zero-centered 
Lorentzian with characteristic frequency 1.75±0.33 Hz. In 
the high-energy band, the power spectrum is noisy but 
compatible with the same shape. 

6.4. U1/U2 power spectrum 

The average power spectrum for the Ul and U2 groups 
appear complex (see Fig. 11211 . Since the exact shape could 
be the result of our averaging of different power spectra, 
we do not present detailed fits here. The Ul power spec- 
trum shows a clear band- limited noise component, with an 
excess around 10 Hz. In the high-energy power spectrum, 
this excess is clearly visible as a broad QPO. The total 
0.1-64 Hz fractional rms is 4.3%. The U2 power spectrum 
is weaker (2.3% rms) and a clear QPO is present around 5 
Hz, in addition to a band-limited noise component. A weak 
excess around 10 Hz is visible. In the high-energy band, 



the shape of the continuum components is similar, but the 
5 Hz peak disappears, while the 10 Hz excess shows as a 
clear QPO peak. 



7. The left branch (bottom section) 

From observation #110 (MJD 52560.41), the PCU2 count 
rate goes below 560 cts/s and only weak power-law noise 
or no noise at all is detected in the power spectrum. From 
Fig. |21 it is evident that all these observations are at rela- 
tively low hardness, and a much smoother path is followed 
mostly monotonically downwards, in the HID. We there- 
fore accumulated one average power spectrum from these 
observations (#110 to #148f). The power spectrum is well 
described by a single power law with slope 0.82±0.02. The 
total integrated rms between 0.1 and 64 Hz is only 1.6%. 
No high-frequency cutoff is observed, with a lower limit 
of 25 Hz. The average power spectrum along this branch 
is similar to that of the weak-signal power spectrum de- 
scribed above, although the power-law here is slightly flat- 
ter. 



8. The bottom branch: back to the start line 

The last observation on the left branch, with no signifi- 
cant signal in he power spectrum is from MJD 52693.73 
(Obs #148f). Just before this date there is a consider- 
able hardening (see Fig. ^ and EI- In the following obser- 
vation (#149, MJD 52694.92), band-limited noise, often 
with a significant QPO, appears once again in the power 
spectrum. The transition is marked with a dotted line in 
Figures Qand Eland corresponds roughly to a hardness of 
0.2. 

The hardening is almost monotonic throughout the 
bottom branch, with the exception of a few observations 
that indicate a short period of softening. During three of 
these softer observations (Obs. #154), the hardness went 
back below 0.2 and no noise components were detected in 
the power spectrum, indicating that the source was back 
in the left-branch region. For all other observations after 
the dotted line, noise was observed in the power spectrum. 
The evolution of the power spectrum can be followed once 
more by examining a plot of integrated fractional rms ver- 
sus hardness as in Fig. Figure El shows this diagram, 
where black symbols indicate observations from the bot- 
tom branch and gray symbols show the same points as in 
Fig.ni 

All observations indicated with circles show flat-top 
noise in the power spectrum, sometimes with a QPO. 
The three observations with hardness below 0.2 show very 
weak noise. One observation, located in the 'type B' re- 
gion from the left branch, shows a very peculiar power 
spectrum consisting essentially of a single peak at 1 Hz. 
These three classes are discussed 
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Fig. 13. Plot of integrated fractional rms versus hardness 
for the observations of the bottom. Open symbols rep- 
resent the same points as in Fig. I1UI right branch and 
type C (circles), diamonds (type B), triangles (U1/U2), 
squares (type a), and dots (weak signal). Filled symbols 
correspond to different power-spectrum shapes observed 
in the bottom branch (see text): flat-top noise (circles), 
weak noise (triangles), 1Hz QPO (diamond). 



Fig. 14. Top panel: 900s of the light curve of Obs. #149h 
(16 second bin size). Bottom panel: corresponding spectro- 
gram (time integration 16 seconds, time step 4 seconds), 
where darker points correspond to higher power. Inset: 
average power spectrum from the full observation. 



8.1. Weak noise power spectra 

These three observations with hardness below 0.2 (Obs. 
#153c, 154, 154b) show very little power so that it is dif- 
ficult to determine the shape of the power spectrum at 
these low count rates. The upper limits on the rms are 
consistent with the value from Sect. 7. 

8.2. 1Hz QPO 

Observation #149h, on MJD 52707.91 appears completely 
different from the surrounding ones. The power spectrum 
and the corresponding spectrogram are shown in Fig. 
1141 The power is dominated by a single QPO peak at 
0.91±0.04 Hz, with a FWHM of 0.33±0.06. The spectro- 
gram shows that this peak is not stable in time and shifts 
from 0.7 to 1.2 Hz over the observation. Aside from the 
frequency, this behavior is similar to that of the type- 
B QPOs described above, including the position in the 
color-hardness diagram. We extracted the QPO phase lag 
as described in the previous sections, but the results is 
consistent with no lags (-0.077±0.076 rad). The fractional 
rms in the QPO peak is ~7%. 

8.3. Flat-top noise power spectra 

All other observations of the bottom branch, which at the 
end bends downwards and becomes vertical, joining the 
initial right branch, show rather strong noise components. 
Some examples are shown in Fig. 1151 

We fitted the power spectra with the same model 
adopted for the top-branch and right-branch observations 



at the beginning of the outburst. Only a maximum of three 
components were needed (see Table 4 and Fig [TB|) . A nar- 
row QPO was found to be significantly present in 10 obser- 
vations. The first observations, with one exception, show 
only one flat-top noise component. When the characteris- 
tic frequencies become low, a second noise component ap- 
pears. Only in two cases it was possible to detect a third 
broad component. Unlike the case of the detections sum- 
marized in Table 2, the QPO frequency is substantially 
higher than the characteristic frequency of the first noise 
component, suggesting that the latter should be identified 
with Lb and not with Lh as in the previous case. This 
is confirmed, although with some deviations, by the plot 
in Fig. 1171 which shows the corresponding WK correla- 
tion. The three noise components L\^Li and L3 in Table 
4 are therefore identified as L&, Lh and Lg. The latter of 
these three identifications can only be tentative, as it is 
not possible to check other correlations for lack of appro- 
priate pairs of frequencies. Comparing Fig. 0] and Fig. 1161 
one can see that an alternative identification could be the 
same as for the top branch. The evolution of v\ as a func- 
tion of count rate and hardness for the bottom branch is 
shown in Fig. [7| (empty circles). Clearly, it is not possible 
to obtain a firm identification for these components. 

The integrated fractional rms also increases with in- 
creasing characteristic frequencies: the Lf, component goes 
from 8% up to 31% at the last observation, while the Lh 
component goes up to 23%, values similar to those at the 
beginning of the outburst. The Lq component has an in- 
tegrated fractional rms of 4-8%. For these observations at 
low count rate, the signal is too weak to measure phase- 
lags significantly different from zero. 
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Table 4. Main characteristic frequencies (in Hz) for the bot- 
tom branch of the HID 



Obs # 


VI 


V2 VI, 


V Q 


Bottom branch 


149 


2 91+0 21 

Zj . ' I 1 1 VJ a ZJ _1_ 




8 59+0 1 8 


149b 
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i tU ) — i — vj . J- vj 
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VJ ■ T: Zj — 1_ U*UU . . . 
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VJ ■ VI Zj — 1 — VJ ■ VJ 1 
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154e 


3 81+0 49 

'J ■ VJ I 1 — V / . 1 t / 




6.12+0.10 
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161 


0.14+0.03 


1.40+0.39 




161b 


0.13+0.04 


0.85+0.21 




162 


0.12+0.03 


2.11+0.58 




162b 


0.23+0.05 


2.92+1.29 
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0.11+0.01 


1.84+0.32 




162c 


0.18+0.03 
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163b 


0.13+0.03 


1.21+0.04 




163c 


0.23+0.06 
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9. High-frequency QPOs 

No significant high-frequency peaks were seen in the sin- 
gle power spectra, either at low or at high energies. Since 
high-frequency QPOs in other systems have been detected 
when the source was in the IMS/VHS state, we accumu- 
lated average high-energy power spectra from all type-B 
and type-A observations and searched for excesses above 
30 Hz. Fits were performed to the power spectrum with 
a model consisting of a flat level (to account for the 
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Fig. 15. Four representative power spectra from the 
bottom-branch observations. The time sequence here is 
bottom to top: observations #161, 158b, 155, 149g. The 
power spectra are plotted in the vP v representation 
(Belloni et al. 1997). 
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Fig. 16. Same as Fig.0J but for the bottom-branch obser- 
vations. 

Poissonian component) plus a Lorentzian. The search was 
done fixing four values of Q: 5, 10, 15 and 20, and with dif- 
ferent starting centroid frequencies in order to make sure 
the best fit was reached. No significant peak was detected. 
The 3cr one-trial upper limits for type A are 2.6%, 2.3%, 
2.0% ad 1.3% for the four values of Q. For type B they 
are 2.6%, 2.3%, 2.0%, 2.8%. 

10. Discussion 

10.1. The outburst 

The results presented above outline a rather clear picture 
for the outburst of GX 339-4. From the combined tim- 
ing and color properties, we can identify four main states, 
which have a correspondence with those determined from 
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0.1 1 10 

v b (Hz) 

Fig. 17. WK correlation (see Wijnands & van der Klis 
1999) for the observations on the bottom branch of the 
HID (filled circles for v\, and empty circles for Vh). The 
small points are from Belloni, Psaltis & van der Klis 
(2002). 

a detailed spectral analysis (Homan et al. in prep.) and 
optical/nIR data (Homan et al. 2005a). The transitions 
between these states correspond to precise locations in 
the HID and can be very fast. In the following we exam- 
ine these states in terms of their timing features. 

10.1.1. The right branch: Low/Hard State (LS). 

The observations on the right vertical branch in Fig. 
12 show hardness and timing properties typical of the 
"canonical" LS. The PSD is dominated by strong band- 
limited noise (see top four power spectra in Fig.|3J). Notice 
the similarity between the second power spectrum from 
the top in Fig.EJand that shown by Belloni et al. (1999) for 
the LS of GX 339-4. This noise can be decomposed into a 
relatively small number of Lorentzian components, which 
have clear counterparts in the power spectra of other sys- 
tems (see Nowak 2000; Belloni, Psaltis & van der Klis 
2002). As usual in black-hole transients, this corresponds 
to the earlier part of the outburst, where the X-ray flux 
increases steadily (see Cui et al. 1999). Interestingly, the 
right branch in the HID is almost vertical, indicating that 
the spectral shape does not change much during the rise 
(see Vignarca et al. 2002; Homan et al., in prep.). In a few 
of the power spectra along this branch, a low-frequency 
QPO is observed (see Grebenev et al. 1999; Belloni et al. 
1999). All characteristic frequencies of the noise compo- 
nents increase as the source brightens. Hard lags increas- 
ing with time are observed along this branch, both for the 
continuum and the QPO (Figs. and El see also Nowak, 
Wilms & Dove 1999). 

The origin of the strong aperiodic variability observed 
in this state is not known. Our results confirm that in 
transient systems the Low/Hard state is associated to the 
early and late stages of the outburst, but not necessarily 



with low luminosities. The observed time evolution of flux 
and characteristic frequencies, both at the beginning and 
at the end of the outburst, suggests that they are cor- 
related with mass accretion rate. The similarity in timing 
and color parameters indicate that the outburst ends with 
similar physical conditions similar to the beginning more 
than one year earlier. 

10.1.2. The top branch: Hard Intermediate State 
(HIMS). 

From observation #25 (MJD 52398), the hardness starts 
to decrease at a faster rate and the source enters the top 
branch in the HID (see Fig. [21). This horizontal path in 
the HID is followed considerably faster in time than the 
LS branch. This softening is associated with an increase 
of the flux from a soft thermal component, with a simulta- 
neous steepening of the hard component (Homan et al. in 
prep.). As mentioned above, the transition from the right 
to the top branch corresponds to a marked change in the 
behavior of the timing parameters (see Fig. 0} and at the 
same time it corresponds to a clear change in the prop- 
erties of the IR/optical/X-ray correlations (Homan et al. 
2005a). 

During this period, a clear type-C QPO appears in the 
power spectrum, which can be linked to that detected in 
some observations of the right branch. The overall noise, 
however, can still be fit with the same components, whose 
characteristic frequencies increase faster with time than on 
the right branch, and whose fractional rms is reduced (see 
Fig.^J. The continuum and QPO phase lags are still pos- 
itive (hard lags soft) and higher than in the LS (see Figs. 
and EJ • This provides a clear link between this branch 
and the LS. This power-spectrum shape has been tradi- 
tionally observed in the Very High State and Intermediate 
State of BHCs (see Belloni et al. 1997; Mendez & van 
der Klis 1997), starting from the early observations of 
GX 339-4 and GS 1124-683 with Ginga (Miyamoto et 
al. 1991,1993,1994), when these timing properties could 
be associated to the presence of a relatively strong con- 
tribution from the hard component, while at softer X-ray 
colors, a different power spectrum was observed (see be- 
low). The power spectrum of the C state of GRS 1915+105 
are clearly associated to this state (see Morgan et al. 1997; 
Trudolyubov 2001: Reig et al. 2000). Notice that along this 
branch, the source does not reach the soft vertical branch 
yet (see Fig. |2J. The end of this branch is again marked 
by a clear transition in the timing parameters. 

We call this branch the Hard Intermediate State (see 
also Homan & Belloni 2005). The fact that the timing 
components appear to be the same as some of those ob- 
served in the Low/Hard State suggests a similar physi- 
cal origin, while the higher values for the characteristic 
frequencies are probably associated to a smaller value of 
the radius of the accretion flow associated to this time 
scale. This indicates that the same spectral component is 
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responsible for this variability (see also Homan et al. in 
prep.). 

10.1.3. The Top/left branch: Soft Intermediate State 
(SIMS). 

From Obs. #34, as the source moves further left on the 
top branch, the timing properties of GX 339-4 change 
abruptly. On Obs. #34, no flat-top noise component is 
detected, and a type-B QPO appears at ~6 Hz (Nespoli 
et al. 2003). It is clear that a rather sharp transition took 
place here. Notice that the difference in hardness between 
Obs. #33 and #34 is small. From this observation on, GX 
339-4 moves irregularly in the HID, remaining in the color 
range 0.05-0.25 (see FigureEJ. No type-C QPO is detected 
along this branch, although we cannot at present exclude 
that the U1/U2 QPOs are of type C. 

This branch is not as straight and monotonic as the 
previous ones, as the source moves up and down in flux 
by a large amount (see Fig.^l, and also its hardness does 
span a rather extended range in log space. Here, after 
grouping of observations together (see Sect. 6), we can 
identify a clear dependence of timing properties on hard- 
ness (see Fig. Above HR=0.27, the power spectra 
show type-C QPOs, in the HR range 0.19-0.27 type-A 
and type-B QPOs, below HR=0.19 are the Ul/U2/no- 
signal power spectra. All these features are observed over 
a large range of count rates, indicating that the spectral 
shape and not source flux is what determines the tim- 
ing properties. Figure HTH indicates that the type A/B/C 
oscillations also correspond to different non-overlapping 
ranges of total 0.1-64 Hz rms. Notice that in the five ob- 
servations where type-B QPOs are observed, the QPO is 
not present all the time. In the intervals without type-B 
QPO, in the observations where the extraction of no-QPO 
power spectra was not problematic, the no-QPO hardness 
is slightly lower and the corresponding integrated frac- 
tional rms drops to ~3%, a value that brings these obser- 
vations in the type-A region of Fig.^J The identification 
of these intervals with type A is strengthened by the de- 
tection of a ~8 Hz QPO in some of those power spectra. 

The characteristics of the type-B QPOs have been ob- 
served before by Ginga in the power spectrum of GS 1124- 
683 and GX 339-4 (see Miyamoto et al. 1991; Takizawa et 
al. 1997). Also fast transitions between two different types 
of power spectrum, corresponding to type-C and type-B 
QPOs, have been observed in those systems (being part of 
the original definition of Very Hard State) , together with 
spectral changes. Fast transitions have been observed with 
RXTE in the bright transient XTE J1859+226 (Casella et 
al. 2004). Both type-B/type- A and type-B/type-C tran- 
sitions were observed: here a very sharp threshold in 
count rate was observed to correspond to the transitions. 
Recently, the transition between the top and the top/left 
branch has been identified as marking the time of the 
launch of relativistic outflows (see Fender, Belloni & Gallo 
2004). All these properties, together with the properties 



of the energy spectra of GX 339-4, lead to the conclusion 
that the two branches correspond to two different physi- 
cal states of accretion. We call this top/left branch Soft 
Intermediate State (see Homan & Belloni 2005). Notice 
that this state, unlike the others, is characterized by a 
collection of different properties in the power spectra 

Finally, softer observations along this branch are more 
difficult to classify. Our Ul/U2/weak subdivision is the 
result of extensive averaging, due to the shortness of most 
observations. 

10.1.4. The Bottom/Left Branch: High/Soft state 

Below a PCU2 count rate of 500 cts/s, corresponding to 
MJD>52559, only weak variability is observed. Looking at 
Fig. |21 one can see that from this time on the movement in 
the HID is monotonically decreasing in count rate, with a 
general softening. We identify this branch with the canon- 
ical High/Soft State. A long RXTE/PCA observation to 
GX 339-4 in this state during a previous transition showed 
a similar power spectrum with a slightly flatter power-law 
slope 0.62 (Belloni et al. 1999). Notice, however, that the 
observations of the upper left branch labeled "weak sig- 
nal" in Fig. (see Sect. 6.3) have properties that are also 
consistent with the High/Soft State 

10.1.5. The bottom branch: back to the 
Hard-Intermediate and Hard states 

The bottom branch corresponds to a final hardening of the 
spectrum. As one can see from Fig. El the hardening starts 
on the last observation of the left branch, but clear timing 
changes are observed only when the hardness goes above 
~0.2. What we called bottom branch is clearly complex. 
The first part of the branch is horizontal: the source moves 
to the right, although it makes a significant excursion back 
to the region HR<0.2. In this region, once more the power 
spectrum shows little noise level in the form of a power 
law, indicating that it was a brief interval back into the 
High/Soft state. In the rest of the branch however, the 
power spectrum shows clear flat-top band-limited noise 
whose characteristic frequencies follow the same correla- 
tions as those from the top and right branch. Also, Fig. 1131 
indicates that these power spectra show indeed the prop- 
erties of the Hard Intermediate and Low/Hard States and 
fill the gaps in FigEB suggesting a possible connection 
with the U1/U2 points (see Fig. 113(1 . The major differ- 
ence is that here the Lb component clearly appears in the 
power spectrum, so that the source follows the main WK 
correlation (see Fig. I17fl . 

Notice that the bottom branch bends smoothly and 
becomes vertical, ending up parallel to the initial right 
branch. However, here we do not find a sharp transition in 
some parameter and do not know where to put the tran- 
sition to the Low/Hard state, which clearly takes place. 
It is interesting however to note that around MJD 52740 
Obs. #157b) two things take place at the same time: the 
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source reaches the same hardness value at which it left the 
LS (see Fig. [I), and the J-band light curve clearly shows a 
change (see Bailyn & Ferrara 2004). We identify this date 
as the marking of the transition to the Low/Hard State 
(see dotted line in Fig. |2J . Finally, our last point is almost 
exactly in the same location on the HID as our first point. 

It it important to notice that, unlike what happened at 
high count rate in the early part of the outburst, in moving 
from the High/Soft state to the Hard Intermediate State, 
there are no observations showing clear timing features 
that would indicate the Soft Intermediate state, such as 
for instance type- A or type-B QPOs. Interestingly, there 
is one observation (Obs. # 149h) which shows completely 
different timing properties: here the power is dominated 
by one single QPO peak at ~1 Hz, which shows rather 
large variations in time (see Fig. 114(1 . Although no sig- 
nificant phase-lag estimate could be obtained, we could 
speculate that this is a variant of the SIMS observed at 
low luminosities. Clearly, more observations of this type 
would be needed in order to confirm this. A detailed anal- 
ysis of this transition in a number of BHC is presented 
by Kalemci et al. (2004), who also conclude that the tim- 
ing properties change more abruptly than spectral (and 
therefore hardness) properties. A 1-Hz QPO with similar 
properties was observed from XTE J1650-500 (Rossi et 
al., in preparation), and is also present in RXTE archival 
observations of GX 339-4 

10.2. Source states 

In summary, the picture that emerges from the 2002/2003 
outburst of GX 339-4 (see Fig. [21 Homan et al. in prep, 
and Homan & Belloni 2005) is consistent with that com- 
ing from a full spectral analysis and from multiwave- 
length correlations (Homan et al. 2005a). Four main re- 
gions arc identified in the HID, corresponding to four sepa- 
rate states. These states and their transitions can be iden- 
tified from the changes in the parameters of the fast time 
variability. The Low/Hard State (the right branch) corre- 
sponds to the first phase of the outburst and shows little 
spectral variations over a rather wide range of PCA count 
rate. Its timing properties are consistent with those usu- 
ally seen in this state in GX 339-4 and in other systems. 

The High/Soft State (the left branch) has a much lower 
hardness and small amplitude fast variability. More inter- 
esting is what happens in between these two branches, 
which in our data clearly corresponds to transitions be- 
tween the two states mentioned above. The evolution of 
timing properties in the Hard Intermediate State (the top 
and bottom branches) can be seen as a continuation of 
that in the Low/Hard State and is rather similar for the 
high-flux and low- flux branches, corresponding to oppo- 
site transitions. A clear marker for the hard-to-soft tran- 
sition can be found in the timing parameters as well as in 
the IR/X-ray correlations, but the situation seems to be 
smoother for the soft-to-hard transition. Notice that this 
transition appears to be reversible, as for some time GX 



339-4 goes back to the High/Soft state before continuing 
to harden to reach the Low/Hard State. During this state, 
type-C QPOs are observed. The phase lags of the type-C 
QPOs observed here are always positive, even when the 
centroid frequency of the QPO approaches 10 Hz. This is 
in contrast to what observed in GRS 1915+105 and other 
transient systems (see Casella et al. 2004), indicating that 
the phase-lag behavior cannot be categorized in a simple 
way (notice that the contribution of the continuum to the 
phase lags can complicate its measurement, and along the 
top branch the Lh component dominates it, unlike in other 
systems) . 

A fourth region (or branch) in the HID can be iden- 
tified (the top/left branch), possibly present only in the 
high- flux part of the HID. Here the timing properties are 
markedly different from those in the Hard Intermediate 
State, indicating a transition to one of the collective prop- 
erties which we call Soft Intermediate State. Type-A and 
type-B QPOs appear and show rapid variations and tran- 
sitions on short time scales. The transition to this state is 
very sharp and, unlike in XTE J1859+226, no more type- 
C QPOs are observed after it (see Casella et al. 2004) until 
the final part of the outburst. Interestingly, the transition 
into this state appears to be associated with the ejection 
of relativistic jets (Gallo et al. 2004; Fender, Belloni & 
Gallo 2004), which again makes it an important state to 
study. Up to now, it has been observed in a number of 
systems: GX 339-4 (Miyamoto et al. 1991,1993; Takizawa 
et al. 1997; Nespoli et al. 2003; this work), GS 1124-683 
(Miyamoto et al. 1994), XTE J1550-564 (Homan et al. 
2001; Remillard et al. 2002), XTE J1859+226 (Casella 
et al. 2004), H1743-322 (Homan et al. 2005b), and 4U 
1543-47 (Park et al. 2004), making it a rule rather than 
an exception. Our work shows that the behavior in this 
state is more complex, and additional unclassified QPOs 
appear (U1/U2, see above). Power spectra with similar 
features have been observed in other systems such as XTE 
J1550-564 (Remillard et al. 1999), 4U 1630-47 (Dieters et 
al. 2000; Trudolyubov, Borozdin & Priedhorsky 2001) and 
possibly GRS 1915+105 (Morgan et al. 1997; Trudolyubov 
2001). Clearly, these features deserve being studied in 
more detail with longer observations. 

Interestingly, a type-B QPO has not been observed up 
to now in GRS 1915+105, despite the wide range of char- 
acteristics shown by this system (see Fender & Belloni 
2004 for a review). Fender, Belloni & Gallo (2004) inter- 
pret the X-ray/radio correlations of GRS 1915+105 and 
other transients in the framework of the same model; how- 
ever, the time scales involved in GRS 1915+105 are much 
shorter and it is possible that the source moves too fast 
through its HID to show the timing features of this state. 
A full comparative analysis of a number of transient sys- 
tems will be presented in a forthcoming paper. 

The state paradigm presented here is not the same as 
that outlined in McClintock & Remillard (2005). Although 
the Low/Hard and High/Soft state can obviously be iden- 
tified with the Hard and Thermal-Dominant states, our 
two additional states are not the same as the Steep-Power- 
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Law and Intermediate-State of McClintock & Rcmillard 
(2005). A precise comparison is not possible, since their 
classification is based on model-dependent parameters 
coming from spectral fits, but the variety of power spec- 
tra shown to be associated to the SPL in McClintock & 
Remillard (2005), which includes both type-A/B and type- 
C QPOs, excludes a one-to-one relation to ours. 

10.3. Quasi-Periodic Oscillations 

With these results, the A/B/C classification of low- 
frequency QPOs in black-hole candidates is strengthened. 
There are still undecided cases such as our U1/U2 power 
density spectra, but they all correspond to low signal-to- 
noise signals that might turn out to belong to one of the 
three types once observed with higher statistics. 

Type-C QPOs are the most common and are observed 
in all systems that show LS and HIMS (see McClintock 
& Remillard 2005; van der Klis 2005). Their properties 
are rather well defined. They are clearly associated to the 
hard component in the energy spectrum (see e.g. Casella 
et al. 2004) and appear during states when the hard com- 
ponent shows a high-energy cutoff and is possibly asso- 
ciated to an outflow (see Fender, Belloni & Gallo 2004). 
Their frequency is strongly correlated with other observed 
properties such as the source flux and X-ray color. They 
are associated to the QPOs observed sporadically in the 
LS, and a similar association can be made between the 
broad-band noise properties of HIMS and LS. There is no 
consensus as to the physical origin of type-C QPOs, but 
it is clear that these oscillations are particularly impor- 
tant for our understanding of the accretion flow, as they 
relate to a spectral component that could be associated 
to an outflow from the system (see Kanbach et al. 2001; 
Rodriguez et al. 2004). 

QPOs of types A and B are detected over a much 
smaller range of frequencies (see also Casella et al. 2004) 
and, while having a similarly hard energy spectrum, are 
most likely associated to a different hard component, for 
which no high-energy cutoff is observed (see Zdziarski 
et al. 2001). Rodriguez et al. (2004), using simultane- 
ous RXTE/Intcgral observation, showed that for GRS 
1915+105, the energy dependence of the type-C QPO indi- 
cates the presence of a high-energy cutoff, while the energy 
spectrum is composed of two separate components. These 
type-A/B oscillations appear in soft intervals of the HID, 
where no compact jet emission is observed in the radio 
(Fender, Belloni & Gallo 2004). The fast transitions ob- 
served between these two types of QPO and in particular 
transitions involving also the type-C QPO (Casella et al. 
2004) could be indicative of short transitions to the HIMS. 
These mini-transitions could be responsible for small jet 
ejections similar to those observed on short (hour) time 
scales in GRS 1915+105 (Fender & Belloni 2004; Fender, 
Belloni & Gallo 2004) . These transitions are indeed related 
to variations in hardness (GX 339-4) or intensity (XTE 
J1859+226). Their limited range of frequencies, especially 



of the type-B QPOs is particularly interesting. Since a 
transient oscillation at a similar frequency has been ob- 
served in the neutron star system 4U 1820-30 (Belloni, 
Parolin & Casella 2004), if this association is correct the 
frequency of type-B QPOs should have a weak dependence 
on the mass of the compact object. Notice that the 1 Hz 
QPO detected near the end of the outburst appears at a 
hardness value of 0.26, comparable to that at which the 
6 Hz type-B QPOs are seen. Since the properties of this 
1 Hz oscillation (besides its centroid frequency) are com- 
patible with those of type-B QPOs, it is possible that the 
small range of the frequency of these QPOs is due to their 
appearing only at relatively bright phases of the outburst. 
In other words, there is also a count rate dependency of 
the frequency of typc-B QPOs, but most of these features 
appear at a high count rate. In our data, there is some ev- 
idence of a count rate dependence of the 6Hz QPO, as it 
is in the case of XTE J1859+226 (Casella et al. 2004). At 
any rate, the fast transitions between QPO types observed 
here (and in other systems) indicate that state transitions 
take place on very short time scales and are clearly traced 
by the timing properties. The three types of QPO de- 
scribed here (A/B/C) are rather different and are possibly 
tracers of separate spectral components. 

11. Summary and conclusions 

We analyzed a large set of RXTE/PCA observations of 
GX 339-4 during its complex 2002/2003 outburst. From 
the timing and color results, we identify four main states 
of the source, with a distinctive pattern of evolution of the 
outburst. The timing characteristics during the outburst 
are complex, but show simultaneous spectral transitions 
that allow a clear categorization. These states (in time 
sequence) are: 

— Low/Hard state: this state is associated with relatively 
low values of the accretion rate, i.e. lower than in the 
other bright states, although it is clearly not limited to 
low flux intervals. The energy spectrum is hard and the 
fast time variability is dominated by a strong (^30% 
fractional rms) band-limited noise. Sometimes, low fre- 
quency QPOs are observed. The characteristic frequen- 
cies detected in the power spectra follow broad-range 
correlations (see Belloni, Psaltis & van der Klis 2002). 

— Hard Intermediate state: in this state, the energy spec- 
trum is softer than in the LS. The power spectra fea- 
ture band-limited noise with characteristic frequency 
higher than the LS and a rather strong 1-8 Hz type- 
C QPO. The frequencies of the main components de- 
tected in the power spectra extend the broad correla- 
tions mentioned for the LS. However, a sharp transi- 
tion from the LS can be identified by looking at the 
IR/X-ray correlation (Homan et al. 2005a). 

— Soft Intermediate state: here the energy spectrum is 
systematically softer than the HIMS. No strong band- 
limited noise is observed, but transient typc-A and 
type-B QPOs appear, the frequency of which spans 
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only a limited range. The transition from the HIMS is 
very sharp and marked by the transition from type-C 
to type-A/B QPOs in the power spectra, despite the 
relatively small changes in hardness. This state is char- 
acterized by a variety of timing properties, of which 
type-A/B QPOs are only a subset. 
High/Soft state: the energy spectrum is very soft. Only 
weak power-law noise is observed in the power spec- 
trum. The onset of this state is identified with the sys- 
tematic non-detection of features in the power spec- 
trum. 




Fig. 18. Same Hardness-Intensity as in Fig. [21 but areas 
marking the regions corresponding to the source states 
described in the text. 



After the HS, a transition back to the HIMS is ob- 
served, marked by the sudden appearance of band-limited 
noise with type-C QPOs. No clear evidence of a SIMS is 
observed in these later stages of the outburst, although 
one observation shows a strong 1-Hz QPO which could be 
of type B. The HIMS evolves into a LS, roughly follow- 
ing an evolution opposite of that of the beginning of the 
outburst. These four states are observed in other systems 
and they are clearly separated by sharp transitions that 
are essential for their identification. Their location and 
boundaries on the HID of GX 339-4 are shown in Fig. fTSI 

The phenomenology of the states, their evolution and 
the transitions between them could constitute a firm start- 
ing point for the development of physical models for the 
outbursts of black-hole transients. 
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Appendix A: Type-B QPOs in detail 

In this Appendix, we describe the detailed analysis of the 
four type-B observations found in addition to that re- 
ported by Nespoli et al. (2003). Notice that the fits to 
average power spectra described below are only indica- 
tive, since it is clear that there are major fast variations 
in the QPO parameters on rather short time scales (see 
also Nespoli et al. 2003). 
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Fig. A.l. Same as Fig. for Obs. #41. 



— Observation #41'- From Fig. IA.1I it is evident that 
the light curve (top panel) is characterized by os- 
cillations between two flux levels. From the corre- 
sponding spectrogram (bottom panel) one can see that 
the QPO is present only during the high-flux inter- 
vals, where it is correlated with count rate. From the 
spectrogram, we produced two average power spec- 
tra by selecting all power spectra corresponding to 
a count rate <3400 cts/s and >3500 cts/s respec- 
tively (3 PCUs were on during this observation). In 
the power spectrum corresponding to the low count- 
rate intervals (see Fig. lA.lfl . no QPO is indeed vis- 
ible: a fit with a zero-centered Lorentzian yields a 
characteristic frequency ^=2.6S±0.32 Hz and a frac- 
tional rms of ~8%. In the power spectrum correspond- 
ing to high count rates, a narrow QPO is visible. A 
Gaussian fit yields 2/ 9 =6.45±0.03 Hz and <7=0.52±0.03 
Hz, with an integrated fractional rms of ~5%. A band- 
limited noise component is also seen: a Lorentzian fit 
gives ^b=5.74±0.44 Hz and a fractional rms of ~6.5%. 
Additional power is present at low frequencies (<1 Hz). 
The integrated 0.1-64 Hz fractional rms is similar be- 
tween the low and high count rate power spectra. 

— Observation #42: here the situation is slightly differ- 
ent. The QPO is only visible in the high-flux intervals, 
as for Obs. #41, but the low-flux intervals look more 
like dips in the light curve as they are relatively few and 
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Fig. A.2. Same as Fig. IO for Obs. #42. 

sparse. Again, we extracted two average power spec- 
tra from the 4-s spectrogram (see Fig. lA.2|) . During the 
dips (rate <4200 cts/s for 4 PCUs), no significant QPO 
is visible: a good fit is obtained with a broad Lorentzian 
component with ^ maa; =4.23±0.39 Hz, plus an excess 
at low frequencies. Outside the dips, the narrow QPO 
appears at an average frequency of i/ g =5.79±0.02 Hz 
(Gaussian fit) with cr=0.39±0.02 Hz and with an in- 
tegrated fractional rms of ~5.6%. A subharmonic and 
a second harmonic components are also evident. The 
continuum noise can be approximated with a broad 
Lorentzian component (^ ma2; =6.62±0.34 Hz), with as 
usual some excess at low frequencies. As for Obs. #41, 
the integrated 0.1-64 Hz fractional rms is similar be- 
tween the low and high count rate power spectra. 




Fig. A.3. Same as Fig. EH for Obs. #43. 



— Observation #43: from Fig. lA.3l the presence of a simi- 
lar narrow QPO peak is evident. The observation con- 
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sists of three separate segments of ^3400 s exposure 
each separated by 2200s-long gaps. The QPO is present 
only in the first ~1550 seconds of the second inter- 
val, after which it disappears on a time scale of a few 
seconds (see Fig. IA.3(1 . The QPO is clearly related to 
additional low-frequency noise in the light curve (see 
Fig. lA.3|) . Accumulating the 128s power spectrum cor- 
responding to the interval when the QPO is present (in 
this case the 4s resolution is not needed) , a clear peak is 
present. We fitted it with a Gaussian with centroid fre- 
quency ^ 9 =5.35±0.01 Hz and er=0.38±0.01 Hz, for an 
integrated fractional rms of 7.7%. The continuum noise 
can be fitted with two zero-centered Lorentzians with 
characteristic frequencies 0.02±0.01 Hz and 7.38±1.16 
Hz and fractional rms 3% and 4.2% respectively. The 
low-frequency component is clearly the one responsi- 
ble for the noise visible in Fig. ESI The average power 
spectrum outside the QPO interval is quite different. It 
can be fitted with two zero-centered Lorentzians plus 
a broad QPO peak (see also Nespoli et al. 2003). The 
QPO has a centroid frequency of i/ 9 =6.71±0.14 Hz 
and a A=3.27±0.46, for an integrated fractional rms 
of 2.2%. The two noise components have characteris- 
tic frequencies 0.05±0.01 Hz and 0.94±0.19 Hz and 
fractional rms 1.6% and 1% respectively. Interestingly, 
there is a very narrow excess over the best-fit model 
at 6.6 Hz, but an additional component is not required 
by the fits and nothing obvious is visible in the spec- 
trogram. Here the type-B QPO power spectrum has a 
0.1-64 Hz integrated fractional rms of 9.3%, while for 
the broad QPO power spectrum it is 3.2%. Since the 
hardness increases only slightly between the two inter- 
vals, notice that the 3.2% value occupies the position 
of the 'type A' group, suggesting that the broad QPO 
is of type A (see below). 



first part of the observation ^850s shows no evidence 
of the QPO, which appears suddenly together with 
additional low- frequency noise (see Fig. IA.4|) . with 
a variable centroid frequency once again correlated 
with the count rate. All these properties are consistent 
with those observed in previous observations, despite 
the fact that this observation takes place almost four 
months after Obs. #43 (see Tab. 1). We produced a 
spectrogram with time resolution 16 seconds and av- 
eraged the power spectra corresponding to the pres- 
ence and absence of QPO. The QPO is fitted with a 
Lorentzian with z/ 9 =5.71±0.02 Hz and a A=0.44±0.03, 
for an integrated fractional rms of 6.5%. It is accom- 
panied by a clear second harmonic, and by a band- 
limited noise component with characteristic frequency 
0.07±0.02 Hz and integrated fractional rms ~6%. In 
the first part of the observation, when the QPO is 
not present, the average power spectrum is rather 
noisy: it can be fitted with a broad (A=4.05±1.15 
Hz) Lorentzian peak with a characteristic frequency 
of 7.41±0.43 Hz and an integrated fractional rms of 
2.1%, plus a band-limited noise with characteristic fre- 
quency 0.10±0.04 Hz and fractional rms 1.4%. As for 
Obs. #43, the integrated 0.1-64 Hz fractional rms is 
much lower when the type-B QPO is not present: with 
a value of 3.0%, it brings the corresponding point in 
the 'type A' group, although no significant QPO is 
present. 







] J 

" M in 


W — 


1 Ft U 
1 





t} H 1» K- » H !4> « n» St* 

T.HT "'if* 



Fig. A.4. Same as Fig.EUfor Obs. #97. 



— Observation #97: much later in the outburst, the sharp 
type-B QPO appears again in this observation. The 



